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ABSTRACT 

Radial velocity (RV) observations of an exoplanet system giving a value of Mt sin(i) condition 
(i.e., give information about) not only the planet's true mass Mt but also the value of sin(i) for 
that system (where i is the orbital inclination angle). Thus the value of sin(i) for a system with any 
particular observed value of Mt sin(i) cannot be assumed to be drawn randomly from a distribution 
corresponding to an isotropic i distribution, i.e., the presumptive prior distribution . Rather, the 
posterior distribution from which it is drawn depends on the intrinsic distribution of Mt for the 
exoplanet population being studied. We give a simple Bayesian derivation of this relationship and 
apply it to several "toy models" for the (currently unknown) intrinsic distribution of M^. The results 
show that the effect can be an important one. For example, even for simple power-law distributions 
of Mt, the median value of sin(i) in an observed RV sample can vary between 0.860 and 0.023 (as 
compared to the 0.866 value for an isotropic i distribution) for indices of the power-law in the range 
between —2 and -1-1, respectively. Over the same range of indicies, the 95% confidence interval on 
Mt varies from 1.002-4.566 (a = —2) to 1.13-94.34 {a = +1) times larger than A/y sin(i) due to 
sin(i) uncertainty alone. More complex, but still simple and plausible, distributions of Mt yield 
more complicated and somewhat unintuitive posterior sin(i) distributions. In particular, if the Mt 
distribution contains any characteristic mass scale Mc, the posterior sin(i) distribution will depend on 
the ratio of Mt sin(j) to M^, often in a non-trivial way. Our qualitative conclusion is that RV studies 
of exoplanets, both individual objects and statistical samples, should regard the sin(i) factor as more 
than a "numerical constant of order unity" with simple and well understood statistical properties. 
We argue that reports of Mt sin(i) determinations should be accompanied by a statement of the 
corresponding confidence bounds on AIt at, say, the 95% level based on an explicitly stated assumed 
form of the true Mt distribution in order to more accurately reflect the mass uncertainties associated 
with RV studies. 
Subject headings: planetary systems, techniques: radial velocities, methods: statistical 



1. INTRODUCTION 

As is well known the observational study of exoplan- 
ets began with, and in large part has been based on, 
radial velocity (RV) data which allow a measurement of 
the planet's orbital parameters plus a value of Mt sin(i), 
where Mt is its true mass and sin(i) is the angle be- 
tween the direction normal to the p lanet's orbita l plane 
and t he observer's sight lin e (see iMarcv et al.l (|2003l 
[2005h : [Cimming et all (l200l : IJohnsonl (|2009D 'l. Indeed. 



the very classification of an unseen stellar companion as 
an exoplanet is normally made based on the value of 
MTsin(i), hereinafter designated as Mg, the observed or 
indicative mass. 

It would, of course, be preferable to determine Mt 
itself and avoid the degeneracy with the largely unin- 
teresting random variable sin(i), and our understanding 
of exoplanet systems has been greatly advanced by the 
relatively few cases in which the observations of tran- 
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sit events allows the t wo parameters t o be measured 
separately ( seetC harbo nneau et al.l (j2007f ) and references 
therein and I WinnI (|2010[ ) and references there in). 

Nevertheless, the Mt sin(i) degeneracy does not seem 
too serious because it appears to be so simple and well 
understood. In particular, it seems extremely safe to as- 
sume that i is randomly and isotropically distributed or, 
in other words, that the orientation of the orbital plane 
of an exoplanet in space is independent of the direction 
from which we observe it. However, this isotropic distri- 
bution only describes the prior distribution of i, not its 
posterior one, i.e., not the relevant distribution after it is 
conditioned by the measurement of an Mt sin(i) value. 

In order to specify this isotropic i prior distribution of 
sin(i), consider a longitudinal strip of a sphere between 
9 = i and 9 = i -\- d9 hi polar coordinates. The strip 
extends around 27r in (j) (the azimuthal angle), and the 
surface area of the strip is just 27rr^ svL\{i)d9. The proba- 
bility of a randomly oriented vector piercing through that 
area is then just its fractional area of the whole surface 
of the sphere 



f, = 27rr2 sin(i)d6'/27rr2 



(1) 



The pdf of the inclination angle (assuming random ori- 
entation) is thus sin(i). In other words, the probability 
distribution function of the inclination angle of exoplanet 
falls into this range i io i + d9 is, just sin(i). 



sm{i)/^Jl 



In order to determine the prior probability of sin(i), 
simply consider 

f,di = f,dj (2) 

where j = sin(i). After some algebraic manipulation, 
it is easy to show that fj (the pdf of sin(i) falling 
into a range of sin(i) — dsin(i) and sin(i) + dsin(i)) is 

Up to this point, the analysis is straightforward. How- 
ever, complications arise at the next step, the derivation 
of the posterior distribution of sin(i), because it depends 
on the prior or true distribution o f Mr- As we do no t 
yet know the Mr di st ribution (s ee 'Marcv et alJ (|2005| ): 
I Jorissen et all ()2001l ): lUdrv fc Santos (2007)), this con- 
sideration is not only an important one in principle but 
might also be in practice. The present paper is primarily 
intended to investigate this issue, the posterior distribu- 
tion of sin(i) given an observation of Mq, in some detail. 

Before presenting a Bayesian analysis in the next sec- 
tions, it may be helpful to note that the issue resembles 
familiar complications in interpreting photometric data 
that are conventio na lly called M a lmqui st -type biases (see 
iMalmauistl (fT920f) : lEddingtonl (fl9ll : iHogg fc Turlle^ 
(jl998f l) in some respects. Namely, even if the measure- 
ment errors are symmetric and unbiased (and, in the 
simple cases most often analyzed, also gaussianly dis- 
tributed... but that is not essential), the true brightness 
of an astronomical object is normally more likely to be 
fainter than its measured brightness than it is to be 
brighter. The well known reason is that there are usually 
a larger number of fainter objects than brighter ones on 
which the (symmetrical) measurement errors may act to 
produce the observed brightness. 

However, the considerations for sin(i) which we investi- 
gate in this paper are not related to measurement errors. 
It would be unchanged even if all of the observations in 
question were perfect and ideal. Neither is it a selection 
bias on sin(i) of the sort that was briefly considered as 
an explanation of exoplanet RV discoveri es in their ear- 
liest days (see lBlackl (|1997[ ): iGravl (|1997[) ). Rather, we 
are considering the unavoidable consequences of the com- 
bination of a physical variable, Mt with an unobservable 
stochastic one, sin(i), when conditioned by a measure- 
ment of their product. This is, of course, a classical issue 
in Bayesian statistics. 

Section [5] defines the basic question addressed by this 
paper and gives a very simple illustrative example of 
why it can be an important issue. Section [3] presents 
a Bayesian derivation of the equations needed to answer 
the question for any given distribution of masses for a 
population of exoplanets, and Section |4] presents the re- 
sults of the analysis obtained by assuming various "toy 
models" for the true exoplanet distribution of masses. 
We then discuss observational selection effects briefly in 
Section [5] and conclude in Section [6] with a discussion of 
its practical implications for RV studies of exoplanets. 

2. ILLUSTRATIVE EXAMPLE 

The question we wish to analyze can be formulated in 
two equivalent but slightly different forms, one describing 
the Mt distribution and one the sin(i) distribution: 

1) What is the probability that Mt is less than X, 
given that RV data yield Mq {= MTsin(i))? The an- 



swer may be written as P{Mt < X\Mo) and depends on 
P{Mt), the intrinsic distribution of exoplanet masses. 

2) What is the probability that sin(i) is less than Z, 
given that RV data yield Mq (= Mt sin(i))? This answer 
may be written as P(sin(i) < Z\Mq) and also depends 
on P{Mt). 

To relate the sin(i) probability distribution and the 
true mass distribution, it is simply: 

P{Mt < X\Mo) = 1 - P(sin(i) < Z\Mo) (3) 

To illustrate the fundamental issue, consider the fol- 
lowing toy model: Suppose that all exoplanets have a 
true mass of either 1.0 AIj or 2.0 Mj where Mj is the 
mass of Jupiter and that there are an equal number of 
exoplanets with each of these masses. If an exoplanet 
is determined to have Mq ~ Mt sin(i) — 0.5 Mj, the 
value of sin(i) is obviously either 0.5 or 0.25 depending 
on whether it is one of the low or high true mass exoplan- 
ets, respectively. Moreover, since a sin(i) value of 0.5 is 
about 2.236 times more likely than one of 0.25 for the 
prior (isotropic i) distribution of sin(i), it follows that 
the posterior distribution of sin(z) for this system con- 
sists of two S functions, one at 0.5 and one at 0.25 with 
the former having an amplitude 2.236 times that of the 
latter. 

Since any intrinsic exoplanet mass distribution could 
be arbitrarily well approximated by a series of 6 func- 
tions, we can conclude that the the intrinsic mass distri- 
bution affects the posterior distribution of the sin(i), for 
any particular observed value of il/g . 

3. BAYESIAN DERIVATION OF THE POSTERIOR 
DISTRIBUTIONS 

Consider a planet at any mass Mt- Given that 
the inclination angle i is randomly (isotropically) dis- 
tributed, we know that (from the previous sections) the 
pdf of sin(i) falling into a range of sin(i) — dsin(i) and 



sin(z) + dsm{i) is sin(i)/i/l — sin (i). The culmulative 
probability of sin(i) < Z will then be 

P(sin(i) < Z)= sm{i)/^Jl-sm'^{i)dsm{i) (4) 



which is simply 

P{sin{i) < Z) = l- cos(arcsin(Z)) (5) 

This is however not surprising, since we know (from 
trigonometric argument) that: 



P{i < x) ~ 1 — cos(a;) 



(6) 



And therefore, we know that the prior probability of 
finding sin(i) less than Z, which is equivalent to the prior 
probability of observed mass Mq given Mt is just 



P(sin(i) < Z) = 1 — cos(arcsin(Z)) 



(7) 



This, of course, is simply the prior pdf of sin(i), derived 
in Section 1. This distribution function is plotted in the 
upper panel of Fig. [1] in cumulative form, we can also 
look at the probability of observed mass Mq given Mt, 





Fig. 1. — The upper panel shows the cumulative prior distribu- 
tion of sin(i). Note the long tail toward small sin(i) values and thus 
the significant probability for a planet to have a true mass substan- 
tially larger than its observed value of My sin(i). The lower panel 
shows the cumulative prior distribution of Mq given Mt = 1- The 
two curves (as expected) are the same. 



which is simply the foUowing: 
P(Mo|A/t) = 



Mo 
Mt 



1 _ (MSL]2 



(8) 



Thus, the culmulative prior probabiUty of observed 
mass Mq given Mt (note that we are here fixing true 
mass Mt and calculating the prior distribution of the 
observed mass Mq, while in the analyses which follows 
we will do just the reverse to obtain posterior distribu- 
tions) and setting Mt = 1 is just: 



P{Mq < X\Mt = 1) = 1 - Vl -X2 



(9) 



As expected, the plot has the same behavior as in 
P(sin(i) < Z), as plotted in the lower panel of Fig.[T] 

Proceeding now to posterior distributions, Bayes' The- 
orem states 



Assigning A ~ Mt and B ~ Mq, we directly obtain 



P(MT|Afo) 



P{Mq\Mt)P{Mt) 
P{Mo) 



(11) 



Equation |8] gives the first term in the numerator of 
Eq. 111! The second term in the numerator is an un- 
known function (which ultimately may be determined 
from observations), but it is possible to consider simple 
toy models, plausible guesses and even theoretical esti- 
mates for P{Mt) and thus explore their consequences. 
Finally, the denominator of Eq. [Tl] can be obtained from 

P{Mq) = / P{MQ\MT)P{MT)dMT (12) 

Therefore 

PiMTlMQ) = PJMqIMtWIt) 
^ ' ' J P(MQ\MT)P{MT)dMT 

which provides the desired posterior distribution 



P(MT|Afo) = 



/^^O 



= P(AfT) 



/ 



i\/o 



(14) 
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■.P{MT)dMT 



for Mq < Mt- 

It is frequently most interesting to consider instead the 
cumulative probability distribution at which Mt < X, 
requiring the integration of the numerator of Eq. [13] up 
to X. This gives 



rX 



P{Mt < X\Mq) 



J P{MT\MQ)P{MT)dMT 
J P{MT\Mo)PiMT)dMT 



(15) 



Since if the observed mass is Mq, then the true mass 
Mt has to be larger than or equal to Mq, since Mq — 
MTsm{i) and sin(z) < 1, therefore the lower integral 
limit is Mq, and the upper mass limit could be as large 
as physically possible for mass of a planet (Mmax)- 

Therefore, we have the following: 
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(16) 
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--P{MT)dMT 



The upper bound on the integral in the denominator 
Mmax is somewhat arbitrary, corresponding to the max- 
imum mass of any planet drawn from the P{Mt) dis- 
tribution. However, the value of Mmax affects only the 
normalization of P{Mt < X\Mq), not its form. 

This formulation in terms of Mt most transparently 
displays the underlying logic of the derivation. However, 
the same approach can equally well give the answer to 
question ^2 above, since the two are equivalent. In par- 
ticular, 



P(sin(i) < Z\A'Iq) = P{Mt > X\Mq) 



(17) 



P{A\B) 



PiB\A)P{A) 



(10) 



And, thus 
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(18) 

Due to our current ignorance of the true P{Mt), we 
cannot evaluate these expressions uniquely for the actual 
observed values of Mg of known exoplanets. It is never- 
theless instructive to do so for various assumed P{Mt) 
distributions. We devote the remainder of the paper pri- 
marily to that exercise. 

4. POSTERIOR DISTRIBUTIONS FOR ASSUMED TRUE 
MASS DISTRIBUTIONS 

In order to investigate the size and character of the sta- 
tistical effect under discussion, we will apply the formula 
derived in the previous section to a series of simple "toy 
models" of P{Mt), some of which might turn out to re- 
assemble reality at least qualitatively. From the previous 
derivation, P{Mt) is assumed to be a normalized prob- 
ability, however, in the following discussion, we will be 
adopting /mt (y) ^^ the absolute true mass distribution 
(with no normalization). This is simply due to the fact 
that both of the derived Bayesian equations (Eq. [1^ and 
Eg. fTS]) are normalizable quantities, so using the absolute 
true mass distribution /mt (y) is not a problem. 

4.1. Power Law Mt Distributions 

Beginning with a particularly simple possibility, we 
now assume the distribution of true masses of exo- 
planets follow a power-law, thus we adopt the form 
fMriy) — ^y"i where both A and a are constants. 
Then we can evaluate the the main integral (hereafter 

<i>(Mo,a)=/' 



Mq 



, —^AM^dMr) of Eq. [T6]for several 

cases (a = l,a = 0, Q! = — l,a 
following: 
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(19) 
For the distribution of sin(i), we can refer to Eq|31 
We use this result to plot a few specific cases in Fig [21 
assuming various values of a. 

Note from Fig. [5] that a = —2 gives a median sin(i) 
value of 0.860 while an equally plausible value of a = — 1 
gives a median sin(z) value of 0.704. If the mass distri- 
bution is an increasing function of mass, the resulting 
median sin(j) value will be reduced quite dramatically; 
for example a = 1 gives a median sin(i) value of 0.02! 

It is equally easy to generate the corresponding 
P{AIt < -'^l^^o) distributions, using Eq. [HI as shown 
in Fig[3l Since the observed mass is set to 1, thus, the 
true mass has to be larger than 1, and as the power- law 
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1,Mo=1,M^3,<=100 



Fig. 2. — The posterior probability function P(sin(j) < Z\Mo) 
assuming a power-law mass function. When we have a large a, 
there are many large planets, thus, we are more likely to have a 
small sin(i) in order to match the observed A4o. When we have a 
small a or a negative a, there are many small planets, so we have 
a higher probability of a larger sin(j) to match the observed Mq. 
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Fig. 3. — The posterior probability function P(Mt < X\Mo) 
assuming power-law mass function, with a = 1,-1,-2, plotted 
with an assumption of Mq = 1, which is in arbitrary unit. 

index increases, (which means a larger number of high 
mass planets in the true mass distribution), the proba- 
bility of finding a planet below X decreases (as seen in 

Fig [31). 

4.2. Power Law Plus a Delta Function Mt 
Distributions 

The distribution of planetary masses in the Solar Sys- 
tem, the highly non-linear and at least partially non- 
gravitational nature of planet for mation as well a s soni e 
specific theoretical models (see iKokubo fc Ida! (I1996D: 
iKokubo et an(|2006[) ; IIda fc Linl(|2008D and lBaraffe et all 
( 20101 ) and references therein) suggest that the P{Mt) 
distribution might contain one or more characteristic 
masses, rather than being an entirely scale free power 
law. In order to investigate the implications of such a 



P{A1t), we consider a toy model in which some of the 
exoplanets are distributed in an a = —2 power law pop- 
ulation while the others all have the same mass Ale. We 
may then again evaluate the expressions of Section II 
directly. 

Thus we have fM^iy) = Ay" + BS{y — Mc), where 
Mc is the critical mass scale of interest. It is convenient 
to introduce the dimensionless parameter rj, defined by 
Mq = rjMc and to set A and B equal. Without loss of 
generality Ale = 1 is adopted {i.e., Mc is defined as the 
unit of mass) for purposes of plotting and giving numer- 
ical values. We can then obtain the following: 



P{Mt < X\AIo) 
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(20) 
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given that Mq < Mc < X < Al^ax • It is easy to see that 
the addition of . "^ will increase the probability that 

AIt is smaller than X. 
Furthermore, if Mq < X < Ale < Almax, we have 



PIMt < X\AIo) 
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(21) 
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V^ 



This makes sense as the critical mass scale is not within 
the boundary that we consider {AIt < X), so the prob- 
ability decreases. 

Finally, if Mq < X < Ad„iax < Ale, then the results 
are similar to the original situation when Jmt (j/) = Ay°' 
except that some of the planets are in the S function part 
of the distribution, thus reducing the relative probability 
of sampling the power-law portion: 
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We can also obtain the distribution of sin(i): 
P(sin(i) < Z|Mo) - 



rj V 



Z2 



2 1 



1 



Ml" 



(22) 



(23) 



if < sin(i) < AIo/Mc < Mq/X. 

If we set Z = ^ and AIq = rjAIc, while M^ = 1, 
then we can plot the following figure Fig |4l It illustrates 
the discontinuity in the probability at the delta function 
(i.e., when rj ^ Z). 

4.3. A Solar System Like Mass Distribution 

Turning now to a more complex but also more phys- 
ically plausible distribution, we analyze the case of ex- 
oplanet masses distributed in a way similar to that of 
Solar System planets. This distribution can be mod- 
eled very roughly as two power-laws separated by a 
gap in mass. One power-law lies at a low mass range 
(the terrestrial planets) while the other lies at a much 
higher mass range (the giant planets). We consider a 
toy model with 2 power-law mass distributions, one ex- 
tending from IMc to 20Mc, while the other power-law 
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Fig. 4. — The upper panel shows the posterior distribution of 
P(sin(i) < Z\r]) for a power-law (with a = —2) and a delta function 
at Mc. We varies Mq (observed mass) with respect to Mc (where 
the Dirac delta function is) by varying rj as Mq = rjMc. One can 
see that there is a significant jump in probability of P(sin(i) < Z\ri) 
whenever Mq = Mc (aka rj = Z). The lower panel shows a similar 
plot, except for P{M < X|Mo = rjMc) as we vary X. 

is for 400Afc to SOOOMc. There are no planets in the 
range between 20Mc and 400A'/c. We also assume the 
two power-laws have the same power index, and also 
same coefficient (i.e. Jmt (y) = Ay'^ in range of IMc to 
20Mc and Jmt {y) = Byl^ in range of 400Mc to SOOOMc 
where A = B and a ~ /3.) We plot the probability 
P(sin(i) < Z\ri) as rj varies (the ratio of the observed 
mass AIq to the critical mass Ale) for a = —2 in Fig[Sl 
Note that the probability P(sin(i) < Z\ri) can saturate 
very near either unity or zero over a substantial range of 
Z values depending on the value of 77. 

5. OBSERVATIONAL SELECTION EFFECTS 

In the preceding analysis we have consistently assumed 
that exoplanets discovered by the RV method uniformly 
(i.e., without bias) sample the distribution of AIt and 
sin(z) values in nature. Obviously, this is unrealistic. In 
reality, both variables (and others) influence the proba- 
bility that a given exoplanet system will be detected in 
an RV survey, and this selection bias in turn affects the 
likely values of both Mt and sin(i). 




Fig. 5. — Both of the above panels show the posterior probabil- 
ity distribution of sin(j) and M for a Solar System like mass dis- 
tribution. The upper panel shows the probability P(sin(j) < Z\r]) 
as we change rj, thus the observed mass Afo goes from O.lMc to 
lO.OOOA/c. The lower panel plots P{M < X\Mo = r}Mc) as we 
change X. Mr_ is set to 1 as usual. 

Happily, this complication does not fundamentally al- 
ter our results because the basic effect discussed in this 
paper is a purely statistical one, independent of any ob- 
servational biases. More specifically, one could conduct 
an exactly parallel analysis in which the true distribu- 
tions of Mt and sin(i), which appear in equations Eg. If 61 
to Eq. [ini are replaced with the biased distributions 
which a particular RV survey samples, if its selection 
function can be determined reasonably accurately. 

A very simple example would be a case in which the 
probability of an RV survey detecting an exoplanet of 
mass Mt is given by some selection function S{Mt), 
independent of sin(i) and other properties of the sys- 
tem. In that case, it suffices to replace P{Mt) with 
S{Mt)P{Mt) everywhere it occurs in the equations and 
proceed as before. 



system will not he drawn from its prior distribution, cor- 
responding to an isotropic distribution of i as is often 
assumed, at least implicitly. 

The relevant, i.e., posterior, probability distribution 
of sin(i) depends sensitively on the distribution of true 
masses Mt and the observed mass A/q — MTsnv{i). 
Since the former is not well constrained, either empir- 
ically or theoretically, at present the true mass Mt of 
such a system cannot be trivially estimated from the 
valu e of Mq as is also oft en assumed to be the case 
seelButler et all (1200 4): M ayor et all (I2005D : iLovis et all 
200(it ): I Wright et al.r(|2008i) : IMavor et al.l (|2OO90 ^. 

This means, for example, that it is difficult to identify 
the least (or most) massive RV exoplanets discovered to 
date because selecting low values of Mt sin(i) from an 
observed exoplanet sample is a way of picking out low 
sin(i) values as well as low Mt values. For some possible 
exoplanet mass distributions the observed objects with 
the lowest observed Mt sin(i) will be dominated by sys- 
tems with small sin(i) values rather than small masses! 

It also implies that the distribution of true exoplanet 
masses is not the same as the distribution of AIq with 
a constant offset in mass, corresponding to the inverse 
of the average or median value o f sin(z) in the sarn - 
ple, as is sometim e s assumed (see iMavor et al.l (|2005[ ): 
iButler et al.l (pOOl i lCmnming et all (|2008f )l. 

The moral of the above analysis is that one should 
be sure to respect sin(i) in RV exoplanet studies. For 
example, we urge that RV observers reporting the value 
of Mt sin(i), typically for a newly discovered planet, also 
report a confidence interval for Mt at some standard 
selected level {e.g., 95%) based on some explicitly stated 
assumption for the true exoplanet mass distribution. 

As a simple illustration, the 95% confidence intervals 
for Mt if Mq = I.O are 1.0017 to 4.566, 1.005 to 27.02, 
1.15 to 85.186 and 1.125 to 94.34 for the simple power-law 
Mt distributions considered in Section 4.1 with assumed 
power- law slopes of a = —2,-1,0 and +1, respectively. 

Although these model dependent upper bounds may 
appear less impressive or exciting than the Mt sin(i) 
value itself, they are a less misleading and thus more 
scientifically informative indication of the actual infor- 
mation on any particular exoplanet 's mass provided by 
RV data alone. 
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6. IMPLICATIONS 

The primary implication of the results presented here 
is that in general the value o/sin(i) for a given exoplanet 
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